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' _j ' Abstract. Dense matter as it can be found in core-collapse supernovae and neutron stars is 

'-pi expected to exhibit different phase transitions which impact the matter composition and the 

' equation of state, with important consequences on the dynamics of core-collapse supernova 

explosion and on the structure of neutron stars. In this paper we will address the specific 
^ phenomenology of two of such transitions, namely the crust-core solid-liquid transition at sub- 

Ch saturation density, and the possible strange transition at super-saturation density in the presence 

of hyperonic degrees of freedom. Concerning the neutron star crust-core phase transition at zero 
I and finite temperature, it will be shown that, as a consequence of the presence of long-range 

^ Coulomb interactions, a clusterized phase is expected which is not accessible in the grand- 

\^ canonical ensemble. A specific quasi-particle model will be introduced and some quantitative 

results relevant for the supernova dynamics will be shown. The opening of hyperonic degrees of 
freedom at higher densities corresponding to the neutron stars core also modifies the equation of 
state. The general characteristics and order of phase transitions in this regime will be analyzed 
in the framework of a self-consistent mean-field approach. 

o 

m 

•'-j 1. Introduction 

rN Triggered by the first data taking of X and gamma satellites, as well as by the improved 

capabilities of radio telescopes and high and ultra- high gamma ray detectors, there has been in 
the last ten years an impressive accumulation of observational data on supernova and neutron 
stars in different environments and at different evolution stages. Among the most important 
recent results one should mention the precise measurement of a very massive two solar-mass 
neutron star [Ij which challenges a number of theoretical equations of state of dense matter, and 
the measurement of the cooling curve of the young pulsar Cassiopeia A [2], which gives strong 
constraints on the super- fluid properties of neutron star cores. 

Understanding these data requires a good modelization of the supernova dynamics and of the 
formation process of neutron stars |3lll], but also a precise and detailed control of the relevant 
microphysics. Sophisticated two and three-dimensional core-collapse supernova simulations show 
that matter in the supernova core explores a very large interval of baryonic densities, ranging 
from about p > 10^° g ■ cm~^ 10~^po to several times the normal nuclear density po , 



temperatures between some hundreds of keV and around 80 MeV, and proton fractions between 
« 0.5 and ^ 0.1. 

Both the equation of state and the neutrino cross sections are strongly influenced by the 
possible presence of phase transitions in these wide thermodynamic conditions. In particular, 
first-order phase transitions lead to a softening of the baryonic pressure, while second order phase 
transitions lead to a dramatic increase of the electroweak interaction probability [5] . This is why 
a control of the phase diagram of stellar matter is of great importance in our understanding of 
these complex astrophysical process. 

In this contribution, we will limit ourselves to the thermodynamic conditions where matter 
can be described by nucleonic or hadronic degrees of freedom and will examine two phase 
transitions which are expected to play a role in supernova physics, namely the crust-core phase 
transition at sub-saturation density and the possible phase transition towards strange matter at 
super-saturation density. 

2. Sub-saturation densities and the quenching of the first-order phase transition 

Since the early days of the equation of state modelization it was clearly recognized that, 
because of the strong similarities between the effective nucleon-nucleon potential and molecular 
interactions, the phase diagram of diluted nuclear matter should present a first-order phase 
transition of the liquid- gas type, terminating at high temperature and density in a critical 
point [6J. This conjecture has been confirmed in all phenomenological [TJ |H] or microscopic [S] 
modelizations of symmetric and asymmetric nuclear matter with realistic effective interactions 
compatible with recent nuclear data. 

In these calculations, nuclear matter is defined as an idealized bulk medium of neutrons and 
protons where the electromagnetic interaction is artificially switched off in order to achieve a 
thermodynamic limit. In the physical situation of stellar matter, however, charge neutrality is 
achieved by the screening effect of electrons on the proton charge, , which induces Coulomb 
effects that drastically modify the liquid-gas phase transition associated to uncharged nuclear 
matter. Because of the very different mass between electrons and protons, the compressibility of 
electron and proton matter is very different. Therefore, in all density and temperature conditions 
relevant for neutron star physics, the electron charge can be safely considered as uniformly 
distributed [10]. This means that any baryonic density fiuctuation (which is correlated to a 
proton density fluctuation because of the repulsive symmetry energy) induces a fluctuation in 
the electric charge. A well known consequence of the resulting Coulomb correlations is that the 
low density phase at zero temperature is not a gas, but a Wigner crystal of nuclei immersed in 
the homogeneous electron background [TT|. At finite temperature, the presence of the electron 
background leads to a quenching of the first order liqui-gas transition, which is replaced by a 
continuous transition through a clusterised medium |12j . This Coulomb frustration effect has 
been confirmed by different microscopic models |13] . A continuous evolution from the Coulomb 
lattice to a homogeneous nuclear fluid, passing through the formation of clusters of different 
sizes and strongly deformed dis-homogeneous structures close to the saturation density has been 
reported, both at zero and at finite temperature. 

These calculations are numerically very heavy and a complete thermodynamic characteriza- 
tion of stellar matter under the frustrated phase transition in a microscopic model has not been 
done yet. Such a task can however be performed in the phenomenological so-called nuclear sta- 
tistical equilibrium (NSE) approaches, which treat the bound states of nucleons as new species 
of quasi-particles [TH d^l dSl IT] . 

2.1. A phenomenological model of dilute stellar matter 

Within NSE, the baryonic component of the stellar matter is regarded as a statistical equilibrium 
of neutrons and protons, the electric charge of the latter being screened by a homogeneous 



electron background. As a first approximation, one can consider that the system of interacting 
nucleons is equivalent to a system of non-interacting clusters, with the nuclear interaction being 
completely exhausted by clusterization. This simple model can describe only diluted matter at 
/9 <C po as it can be found in the outer crust of neutron stars, while nuclear interaction among 
nucleons and clusters has to be included for applications at higher density, when the average inter- 
particle distance becomes comparable to the range of the force. A simple possibility \18\ [19] 
is to take into account interactions among composite clusters in the simplified form of a hard- 
sphere excluded volume, and interactions among nucleons in the self-consistent Hartree-Fock 
approximation with a phenomenological realistic energy functional fitted on a large number of 
measured mass and radii. In this version of the model, the thermodynamics can be calculated 
in the ensemble {I3,fij,p) defined by inverse temperature /? = T~^, isovector chemical potential 
1^1 = tJ'n — fJ'p and baryonic density p = Pn + Pp- The total entropy density is given by 

(P) = 1^ ^r^' (Pf) + ^ 1" ^^"S i^Pci) (1) 

where the density repartition between the clustered pd and unbound pf component p = Pd + Pf 
is uniquely defined by the condition of having a single chemical potential p for both components. 
The unbound nucleons partition sum is calculated in the mean-field approximation: 

In Zf,=] (p) = In Zl^^^^ - /5 In Zl^^^^ + Fe^^) - /SppV, (2) 

In this equation e^^ is the Hartree-Fock energy and the non-interacting part of the partition 
sum can be expressed as a functional of the kinetic energy density Tq for neutrons (q) or protons 
{q=p) : 

q,P 

where Tq are self-consistently derived as a function of the associate chemical potential pq and 
effective mass m*. 

The partition sum of clusters for a total particle number A = Vpd is calculated via a recursive 
relation: 

1 ^ 

(^) = ]4 E ('^aZf5,^r {A -a). (4) 

a=2 

where the weight of the different cluster size is given by: 

uJa = VF J2 exp(-/3/f .)exp(/3i///). (5) 

Here, f^^ is a phenomenological free energy functional including the screening effect of electrons 
in the Wigner-Seitz approximation, and Vp = V{1 — Pd/ Po) is the volume fraction available to 
the clusters ^ [l9] . 

The thermodynamics of the extended NSE model was studied in ref. [19]. As it is well- 
known when studying phase transitions with long range interactions j20], the different statistical 
ensembles are not equivalent and it is very important to use the modelization adapted to the 
physical situation under study. In the case of compact stars, the density profile is imposed by 
the gravitational field and the correct ensemble is the canonical one. 

The phase transition quenching and the effect of the statistical ensemble is illustrated in the 
left part of Figure [T| which shows the cluster distributions at a given temperature. As a function 
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Figure 1. Left side: Cluster number densities as a function of the cluster size a at different 
densities (expressed in g/cm?) of the extended NSE model at T = 1.6 MeV and asymmetry 
chemical potential /i/ = 1.68 MeV. Right side: Comparison between canonical (full line) and 
grand-canonical (dashed line) predictions for the cluster distribution in a specific thermodynamic 
condition relevant for supernova dynamics. Figure taken from ref. [T9|. 



of the density, the average cluster size is continuously changing and can never be assimilated 
to a portion of the fluid phase. The right part of the same figure compares the grandcanonical 
and canonical cluster distribution in a thermodynamic situation relevant for supernova physics. 
We can see that accounting for the correct thermodynamics has important consequences on the 
matter composition. Since the cluster abundances determine the electron capture rate, which in 
turn is a capital ingredient for the size of the homologous core and the cooling process, it is clear 
that a control of the phase diagram is important for supernova physics. It is also interesting 
to remark that the most probable abundances concern not so neutron-rich nuclei which are 
potentially accessible in laboratory experiments: a detailed knowledge of the mass, level density 
and electron capture cross section for these nuclei is thus needed to have a reliable description 
of the equation of state. 

2.2. In-medium effects 

In the extended NSE model presented in the previous section clusters are treated as an ideal 
gas characterized by an empirical energy functional which is modified with respect to nuclei in 
the vacuum only through the excluded volume mechanism. This is an improvement compared 
to standard NSE [U [151 El HZ] iii which the Coulomb screening by the electron background 
is the only interaction accounted for, but still represents a crude approximation of the effect 
of the medium [21]. Microscopic modifications of the cluster binding energies due to Pauli 
blocking shifts have been computed ^ and shown to correctly describe the Mott dissolution of 
deuterons in dense matter, but this heavy numerical approach cannot be extended to massive 
clusters. In the original Lattimer-Swesty model p2] the assumption is made that the effect 
of the medium consists in a modification of the cluster surface tension, the bulk energy being 
unaffected. However some interpretations of multifragmentation data suggest that bulk terms 
might also be affected p3j. 

To microscopically account for in-medium modifications of the cluster binding energies, we 
have realized a systematic set of Hartree-Fock calculations in the Wigner-Seitz cell for various 
isotopic chains with a total neutron number ranging from the stability line up to ~ 3000. The 
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Figure 2. Left part: isoscalar and isovector bulk and isoscalar surface term of the in- medium 
modified mass formula extracted from Hartree-Fock calculations in the Wigner-Seitz cell. Right 
part: comparison between the microscopic Hartree-Fock volume energy per particle and surface 
energy per A^^^ and the parameterization eq.(^6| for three isotopic chains as a function of the 



bulk proton- neutron asymmetry. The model surface term is given by aJg{p) — a'^gy^{p,5)6Q 
see eq. Q. 



resulting microscopic nuclear energy, consisting of a volume and a surface term, is parameterized 
with an in-medium modified mass formula 

Ehf = e"^{pn,Pp){Vws-Vci)- 

where Vws is the volume of the calculation cell and pn,Pp,p,S is the neutron, proton, total 
density and asymmetry 6 = {pn — Pp)/ p of the free nucleons constituting the dripping gas. The 
cluster mass and charge A, Z and the gas densities are extracted through a Wood-Saxon fit 
of the local density. The bulk density pQ and asymmetry 5q = {pon — Pop)/po in the model 
are approximated by analytical functions of A,Z, while Vd = A/pQ. The resulting in-medium 
parameters obtained from a preliminary fit of the microscopic Hartree-Fock results through 
eq. Q are shown in the left part of Figure [2| while the quality of the fit is displayed in the right 
part of the same figure. 

To obtain the fit shown in Figj2]we have employed the ansatz 

af{p) = as(l-£j'' ; aZymiP,^) = ^Ssym (l - (7) 

with 7i = 72 = 2 and with assym = —32 MeV. A dependence of a^gym °^ POi^o is implicit. We 
can see that the in- medium modified liquid drop expression eq. ([g]) provides a reasonably accurate 
description of the energy modification due to the external gas. In agreement with the conjecture 
of ref. |22j, the bulk terms are completely unaffected by the external medium (a^ (p) = ay=cst., 
^VsymiPi^) = avs?/m=cst.), while important effects on the surface properties of nuclei are seen. 
The implementation of such an improved energy functional in the extended NSE model is in 
progress. 
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3. Supra-saturation densities and the strangeness phase transition 

If it is well admitted that hyperonic and deconfined quark matter could exist in the inner core 
of neutron stars, a complete understanding of the composition and equation of state of dense 



matter is far from being achieved. Concerning hyperons, simple energetic considerations suggest 
that they should be present at high density [21] • However, in the standard picture the opening of 
hyperon degrees of freedom leads to a considerable softening of the equation of state [211 [23 [2H] , 
which in turns leads to maximum neutron star masses smaller than the observed masses of many 
neutron stars, in particular the high values obtained in recent observations [1]. This puzzling 
situation implies that the hyperon-hyperon and hyperon-nucleon couplings must be much more 
repulsive at high density than presently assumed [23 [281 [29l [301 [311 [32] , and/or that something 
is missing in the present modelization. 

The generic presence of attractive and repulsive couplings suggests the existence, in a model- 
independent manner, of a phase transition involving strangeness. Let us consider a simplified 
situation where matter is uniquely composed of neutrons and A-hyperons. The energy density 
e = efcin + ^pot at zero temperature is obtained with the energy density functional proposed by 
Balberg and Gal 



epot{Pn,PA) = ^ {aNN + bNN) pfi+CNNpi^'^ +^ IAA/OA + CAApX^ 

+ aAnPnPA + CAn ^".^"^ [pZ + pVi > (8) 
Pn + PA 

Due to the attractive part of the nA coupling, a mixture of the two particle species admits a 
bound state at finite density. This rather general feature of the energy density is found within 
other parameterizations and other models, e.g. the parameterization of the energy density 
functional from G-matrix calculations by Vidana et al. [M]. The existence of a minimum in the 
energy functional for symmetric matter means that such a state represents the stable matter 
phase at zero temperature. On the other side, the vanishing density gas phase is always the 
stable phase in the limit of infinite temperature. This implies that a dilute-to-dense phase change 
in strange matter has to be expected on very general grounds. It is intertesting to remark that 
the minimum, for a certain interval of neutron densities, is also observed as a function of the 
A density. This is again due to the interplay between the attractive A — A interaction required 
by hypernuclear data, and the repulsion at high density required by the two-solar mass neutron 
star observation. This information means that the expected phase change is expected to have 
strangeness contributing to the order parameter. 

The presence of a minimum is however not sufficient to discriminate between a smooth cross- 
over and a phase transition. In the next section we will therefore demonstrate the existence of 
a phase transition by explicitly calculating the nA phase diagram. 



3.1. Neutron-A phase diagram 

First-order transitions are signaled by an instability or concavity anomaly in the mean-field 
thermodynamic potential, which has to be cured by means of the Gibbs phase equilibrium 
construction at the thermodynamic limit. At zero temperature, one thermodynamic potential 
is given by the total energy 

Etot (Pn, pa) = epot (Pn, pa) + Cfcrn (Pn, Pa) + (Pn"in + PA"1a) (9) 

The spinodal region is then recognized as the locus of negative curvature of the energy surface, 
Cmin < 0. The corresponding eigenvector defines a direction in the density space given by 

Pn _ CnA _ Cmin ~ Cw (10) 

Pa Cmin Cfin C\n 

with dj = d'^etot/ dpidpj. This instability direction physically represents the chemical 
composition of density fluctuations which are spontaneously and exponentially amplified in 
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Figure 3. Left part: {n, A) mixture at T = 0: Borders of the phase coexistence region (red thick 
hnes), constant-/x„ paths after Maxweh construction (blue thin lines) and jj-n = fJ-A trajectory 
after Gibbs construction (dark green dotted line). Right part: borders of the phase coexistence 
region of the nA mixture for different temperatures, r=0, 5, 10, 20 MeV. The full circles mark 
the critical points. Figures taken from ref. |35| . 

the unstable region in order to achieve phase separation, and gives the order parameter of the 
associated phase transition. 

The zero temperature instability region of the nA mixture is shown in the left part of Fig. |3| 
We can see that a large portion of the phase diagram is concerned by the instability. We 
can qualitatively distinguish three regions characterized by different order parameters. Below 
nuclear saturation density, we observe an isoscalar pn w pA instability, very close to ordinary 
nuclear liquid-gas, with A's playing the role of protons. Increasing neutron density the phase 
separation progressively changes towards the strangeness ps = —pA direction: the two stable 
phases connected by the instability have close baryon densities but a very different fraction of A's. 
For high p\ we observe the same behavior with the roles of neutrons and A-hyperons exchanged. 
The part of the phase diagram at high neutron density comprises the region physically explored 
by supernova and neutron star matter, which is characterized by chemical equilibrium for 
reactions implying strangeness, ps = l^-n — l-^A = 0. The strangeness equilibrium trajectory 
is represented by a dashed line in Figj3| This physically corresponds to the sudden opening of 
strangeness, observed in many modelizations of neutron star matter (see e.g. [33', '26]). 

Let us stress that the existence of a strangeness phase transition at high density is not a 
general model-independent feature, although many models show it. There are others, for instance 
the G-matrix models of ref. [3l], which do not show an instability in this region. The reason is 
essentially due to the absence of a A A interaction in these microscopic calculations. Owing to 
the lack of reliable information on the hyperonic interactions, which would discriminate between 
different models, we cannot affirm the existence of the strangeness phase transition, related to 
this instability, but, turning the argument around, the presence of this phase transition in a 




Figure 4. Left: A fraction as a function of tlie baryon density witliin the condition ^„ = for 
a nA and, respectively, a symmetric pn = Pp npA mixture at r=0. The full lines correspond to 
the stable mean-field results; the dashed lines illustrate the Gibbs constructions; the horizontal 
lines indicate the relative amount of A-hyperons at the frontiers of the phase-coexistence domain. 
Right: A density for the two coexisting phases as a function of temperature for ps = 0. Figures 
taken from ref. 1351 . 



physical system would allow to learn much about the shape of the interaction. 

Up to now we have only presented results at zero temperature. The extension to finite 
temperature, as it is needed for supernova matter, is relatively simple in the mean-field 
approximation. The appropriate thermodynamic potential at non-zero temperature is given 
by the Helmholtz free energy 

/t (Pn, pa) = etot {Pn, pa) - Ts {pn, Pa) (H) 

where s is the mean-field entropy density. Chemical potentials can be obtained by differentiating 
this expression. Then the Gibbs construction is performed from the combined analysis of PniPn) 
at constant /iA, and pa{pa) at constant pn- 

The phase diagram as a function of temperature is presented in the right panel of Fig. |3j 
This phase diagram exhibits different interesting features. We can see that the three regions that 
we have tentatively defined at zero temperature appear as distinct phase transitions at finite 
temperature. The first phase transition (zone II in Figjs]- right panel) separates a low density 
gas phase from a high density more symmetric liquid phase, very similar to ordinary liquid-gas. 
The second one (zone III in Figjs] - right panel) reflects the instability of dense strange matter 
towards the appearance of neutrons and has an almost symmetric counterpart (zone I in Figjs] 
- right panel) in the instability of dense neutron matter towards the formation of A-hyperons. 
Up to a certain temperature, this latter phase transition is explored by the pn = PA trajectory, 
meaning that it is expected to occur in neutron stars and supernova matter. At variance with 
other known phase transitions in nuclear matter, this transition exists at any temperature and 
is not limited in density; it is always associated with a critical point, which moves towards 
high density as the temperature increases. This means that criticality should be observed in 
hot supernova matter, at a temperature which is estimated as Tc = 14.8 MeV in the present 
schematic model. 

The consequences of these findings for the composition of neutron star matter are drawn in 
Fig. [4| The left panel shows the A fraction 1a = PA/{Pn + Pa) as a function of the baryon 
density under the condition /i„ = p\. The crossing of the mixed-phase region with increasing 
neutron density implies that, as soon as the lower density transition border is crossed, the system 



has to be viewed as a dis-homogeneous mixture of macroscopic regions composed essentially of 
neutrons, with other macroscopic regions with around 25% of hyperons. The extension of the 
A-rich zone increases with density until the system exits the coexistence zone, and becomes 
homogeneous again. The density values associated to the phase transition slightly change with 
the energy functional assumed. To give a simple example, we present on the same Fig|4] a 
calculation where A's are put in equilibrium with symmetric pn = Pp non-strange matter, that 
is eq.Q is still used, but the symmetry energy in the nucleon sector bNN is put to zero [33]. 
We can see that only small differences are observed in the strange phase transition considering 
pure neutron matter or symmetric matter. 

On the right panel of Fig. |4] the A-densities for the two coexisting phases are displayed as a 
function of temperature, again under the condition of /x^ = /xaj physically relevant for neutron 
star and supernova matter. The extension of the coexistence region decreases with increasing 
temperature. The latter finally disappears at Tc = 14.8 MeV. This well illustrates the fact 
that the critical point of the strangeness phase transition moves to higher density, crossing the 
physical line p,n = p\ at a given temperature. 

4. Conclusions 

In this contribution we have analyzed the phase diagram of stellar matter as it is produced in 
core-collapse supernovae and proto-neutron stars in the thermodynamical conditions where such 
matter can be described by nucleonic and hadronic degrees of freedom. 

At baryonic densities lower than nuclear saturation, the first-order phase transition of neutral 
nuclear matter is quenched in stellar matter due to the effect of Coulomb frustration. As a 
consequence, matter is expected to be clusterized even at finite temperature in the whole sub- 
saturation regime. 

The composition and thermodynamics of such clusterized matter has been explicitly worked 
out in the framework of a quasi-particle model where inter-particle interactions are modeled 
through excluded volume. Coulomb correlations with the neutralizing electron gas are 
accounted for in the Wigner-Seitz approximation and residual nuclear correlations are included 
in an effective way through an in-medium modification of the cluster energy functional with 
parameters fitted from Hartree-Fock calculations with a realistic effective interaction. The same 
formalism is used to evaluate the free nucleon self-energies. 

The expectation of the phase transition quenching is confirmed by this model, which shows 
that the single nucleus approximation widely used is a poor approximation at finite temperature 
and a large distribution of cluster species has to be included in supernova modeling. 

In the supersaturation regime, we have used the same mean field approximation with effective 
phenomenological energy functionals in order to describe the characteristics of the phase diagram 
in the presence of strangeness. 

We have shown that already a simple system composed of neutrons and A hyperons in thermal 
equilibrium presents a complex phase diagram with first and second order phase transitions. 
Some of these phase transitions are probably never explored in physical systems. However, 
a possible phase transition at super-saturation baryon densities, from non-strange to strange 
matter is expected to be observed both in the inner core of neutron star and in the dense regions 
of core-collapse supernova. For this latter phenomenology, a critical point is predicted and the 
associated critical opalescence could have an impact on supernova dynamics |1] . The immediate 
consequence of that is that the opening of hyperon channels at high density should not be viewed 
as a continuous (though abrupt) increase of strangeness in the matter, observed in many models 
of hyperonic matter, cf e.g. [26j , but rather as the coexistence of hyperon-poor and hyperon-rich 
macroscopic domains. 

An extension of this simple model to include all possible hyperons and resonances is in 
progress. 
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